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Abstract 

By modehng a dielectric medium with two independent reservoirs, i.e., 
electric and magnetic reservoirs, the electromagnetic field is quantized 
in a linear dielectric medium consistently. A Hamiltonian is proposed 
from which using the Heisenberg equations, not only the Maxwell 
equations but also the structural equations can be obtained. Using 
the Laplace transformation, the wave equation for the electromag- 
netic vector potential is solved in the case of a homogeneous dielectric 
medium. Some examples are considered showing the applicability of 
the model to both absorptive and nonabsorptive dielectrics. 



1 Introduction 



In a homogeneous and nondispersive medium, the photon is associated with 
only the transverse part of the electromagnetic field, in contrast in an inhomo- 
geneous nondispersive medium, the transverse and the longitudinal degrees 
of freedom are coupled. In this case the quantization of the electromagnetic 
field can be accomplished by employing a generalized gauge [1,2], that is, 
V ■ {e{f)A) = 0, where £{r) is the space dependent dielectric function. Using 
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the gauge ^{eij{r)Aj) — [3], the generahzation of this quantization 
to the case of an anisotropic nondispcrsivc medium is straightforward. 
The quantization in a dispersive and absorptive dielectric, represents one of 
the most interesting problems in quantum optics, because it gives a rigorous 
test of our understanding of the interaction of hght with matter. The dis- 
sipative nature of a medium is an immediate consequence of its disspersive 
character and vice versa according to the Kramers- Kronig relations. This 
means that the validities of the electromagnetic quantization in nondissipa- 
tive but dispersive media is restricted to some range of frequencies for which 
the imaginary part of the dielectric function is negligible, otherwise there will 
be inconsistencies. 

In the scheme of Lenac [4], for dispersive and nonabsorptive dielectric media, 
by starting with fundamental equations of motion , the medium is described 
by a dielectric function s{f,u), without any restriction on its spatial behav- 
ior. In this scheme it is assumed there is no losses in the system, so the 
dielectric function is real for the whole space. The procedure is based on an 
expansion of the total field in terms of the coTiplcd cigcnmodes, orthogonality 
relations are derived and equal-time commutation relations are discussed. 
Huttner and Barnett [5] have presented a canonical quantization for elec- 
tromagnetic field inside the dispersive and absorptive dieletrics based on a 
microscopic model in which the medium is represented by a collection of inter- 
acting matter fields and the absorptive character of the medium is described 
by interaction of the matter fields with a reservoir containing a continuum 
of the Klein-Gordon fields. In their model, eigen-operators for coupled sys- 
tems are calculated and electromagnetic field has been expressed in terms of 
them, the dielectric function is derived and it is shown to satisfy the Kramers- 
Kronig relations. 

Gruner and Welsch [6], presented a quantization method of the radiation 
field inside a dispersive and absorptive linear dielectric starting from the 
phenomenological Maxwell equations, where the properties of the dielectric 
are described by a permitivity consistent with the Kramers- Kronig relations, 
an expansion of the field operators is performed that is based on the Green 
function of the classical Maxwell equations which preserves the equal-time 
canonical field commutation relations. 

Suttorp and Wubs [8] in the framework of a damped polarization model have 
quantized the electromagnetic field in an absorptive medium with spatial de- 
pendendence of its parameters. They have solved the equations of motion 
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of the dielectric polarization and the eletromagnetic field by means of the 
Laplace transformation for both positive and negative times. The operators 
that diagonalize the Hamiltonian are found as linear combinations of canoni- 
cal variables with coefficients depending on the electric susceptibility and the 
dielectric Green function. Also the time dependence of the electromagnetic 
field and the dielectric polarization are determined. 

The macroscopic description of a quantum damped harmonic oscillator is 
presented in terms of Langevin equation [9,10]. The coupling with the heath 
bath in the microscopic theory corresponds to two terms in the macroscopic 
equation of motion of a damped harmonic oscillator, a random force and a 
memory function as 

X + / dt'i^U - t')xit') - ul / dt'vit - t')xit') = FnH), (1) 
Jo Jo 

where F/v(t) is the random force and uq is the frequency of the oscilla- 
tor. Matloob [11] quantized the macroscopic electromagnetic field in a linear 
isotropic permeable dielectric medium by quantizing the Langevin equation 
and associating a damped quantum harmonic oscillator with each mode of 
the radiation field. There are some another approaches for quantizing the 
electromagnetic field, for example see [12-22]. 
The structure of this article is as follows: 

In section 2, we quantize the electromagnetic field in a linear dielectric 
medium by modeling the dielectric medium with two independant reservoirs, 
i.e., magnetic and electric reservoirs respectively, where the electric reservoir 
interacts with the displacement vector field through a minimal coupling term 
and the magnetic reservoir interacts with the magnetic field. It is shown that 
both the Maxwell and structural equations can be obtained from Heisenberg 
equations. In section 3, we solve the equations of motion for a homogeneous 
dielectric medium using the Laplace transformation. Finally we give some 
examples indicating the applicability of the model to the case of absorptive 
or nonabsorptive dielectrics. 

2 Quantum dynamics 

Quantum electrodynamics in a linear dielectric medium can be investigated 
by introducing two reservoirs that interact with the electromagnetic field 
through a new kind of minimal coupling. The Heisenberg equations for 
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eletromagnetic field and the reservoirs give not only the Maxwell equations 
but also the structural equations which relate the electric and magnetic fields 
to the electric and magnetic polarization densities respectively. The vector 
potential of the electromagnetic field can be written in terms of plane waves 



A{r,t)= d'kY^ 

A=l 

where oj^ = c\k\ , Eq is the permitivity of the vacuum and e^;^, (A = 1,2) 
are polarization vectors which satisfy 

e^kx-k = 0, (3) 

These recent relations guarantee that the vector potential in Q satisfies the 
Coulomb gauge V ■ A{f, t) = 0. 

Operators aj^xit) '^\x^^^ annihilation and creation operators of the 
electromagnetic field and satisfy the following equal time commutation rules 

[a^,,it),al,^Xt)] = 5ik-k')5xx'. (4) 

The conjugate canonical momentum density of the electromagnetic field 
7ip{f,t) and also the displacement vector operator D{f,t), are by definition 

(5) 

From this definition it is obvious that V ■ D = which is the Gauss law 
in the absence of external charges. The commutation relations (jH lead to 
the commutation relations between the components of the vector potential 
A and the displacement vector operator D 

[Ai{f, t), -D,{r', t)] = [Ai{f, t), 7T,{r', t)] = thS^^{f- P), (6) 

where 6j^{f—r') = J (Pke^'''^^~'^'\6ij — j^), is the transverse delta func- 
tion with the following properties. 

1. Let the transverse F^(f,t) and the longitudinal F^^{f,t) coponents of an 
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arbitrary vector field F(r, t) be defined as 

F^{f, t) = F{f, t)- j d^r'V ■ F{f', t)VG{f, f'), (7) 
Fll(f,t) = J £r'V ■ F{r',t)VG{r,f'), (8) 

where G{f,r'), is tlie Green function 

G{f,?)^-—^, (9) 
47r|r — r\ 

then as the first property, one can easily show that 

KHr,t) = ^ / d'r'6lj{f- r')F^{r'). (10) 

2. The second property is 

T.ir^W-^')-^. i = l>2,3. (11) 
1=1 ^-^i 

The Hamiltonian of the electromagnetic field inside the dielectric can be 
defined as 

J 2eo 2/10 

= j:Jd'kho^,[al{t)a^,,{t) + l], (12) 



A=l' 



where //q is the magnetic permitivity of the vacuum. 

Now we propose the total Hamiltonian, i.e., electromagnetic field plus reser- 
voirs like this 

H{t) = d'r V ^ ' ^^ + ^ ' ^ ' ^ -Vx A f,t .M f,t } 

J 2eo 2yUo 

+ +H, + H^, (13) 

where and H^n are the Hamiltonians of the electric and magnetic reservoirs 
respectively 

He{t) = E / d'q [ d^khw^[dl{k, q, t)d,{k, q, t) + h (14) 



Hm{t) = E / d'q I d'khuMik, q, t)K{k, q, t) + h (15) 

The annihilation and creation operators d,y{k,q,t), dl{k,q,t), b^{k,q,t) and 
bl{k,q,t) of the electric and magnetic reservoirs respectively, satisfy the fol- 
lowing equal time commutation relations 

[d^{k,q,t),dl,{k',q',t)] = 5uu'5{k - k')6{q- 

[K{k, q, t), bl{k', q', t)] = 5,,,5{k - k')5{q- (f). (16) 

Accordingly the electric and magnetic polarization density operators P{f, t) 
and M(r, t), in (fT^. can be expanded as 

P{f,t) = j2 f -4= I d'k[f{u;^,,f^d,{k,q,t)e'^^-^^+r{u^,,f^dl{l 

(17) 

M(r,t) / ^££= / d%g{u^,,^h,{Klt)e'^-' - 9\u:^k^^hl{k 

•—1 ■ I [2tiY 



u=l ' 



(1^ 



where 

vM = ^=|> 

Su{q) = qy. e^g, y = 1,2, 

sM = q- (19) 

In definition of polarization densities (fTTj) . the functions /(cu^, r) and g{ujj^, r), 
are the coupling functions between the electromagnetic field and the electric 
and magnetic reservoirs respectively and are dependent on position r for 
inhomogeneous dielectrics. Using (jU)) and (jlUj) . we can obtain the Heisenberg 
equations for A and D 

?^.l.^^J^r,t)].-£m^pM, (20) 
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dD(f,t) i,~ -^^^ V xV X A(f,t) ^ ,-.|,^ , 
ot n fio 

where and M"*" are transverse components of P and M. Now define 
transverse electrical field E-^, magnetic field B and magnetic induction field 
if as 

i^ = - — , B = VxA, H= M, (22) 

dt fiQ 

respectively, then ()20|) and ()2ip can be rewritten as 

D = eoE^ + P^, (23) 

^ = VxH\ (24) 

in the absence of external charge density, we have = EqE^^ + P" = and 
we can define the longitudinal component of the electrical field as = — 
Combination of and (|7T|) leads to 

1 d'^A dP-^ 

Using commutation relations we find the Heisenberg equations for op- 
erators dfiu{k,t) and bfi^{k,t), 

dyik, q, t) = ^[H, dy(k, q, t)] = 



-iu^,d,ik, q, t) + ^ Z— J d'r'riu^,, r')E{P, t)e-''^-'' ■ v,{q), 



(26) 



— ♦ z — — * 

bu{k,q,t) = -[H,K{k,q,t)] 



luMk, Q, t) + / / d;'r'g*{u,^, r')e-'^'-'' B{r' , t) ■ sM, (27) 
hJ (27r)? 
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with the following formal solutions 

d,{k,q,t) = 4(A?,g,0)e-*"£* + 

rft'e-*^£(*-*') / d^r' f*{ujj:,r')e-'"'-'^' E{r' ,t') ■ v^{q), 



hJ{2n 



|3 Jo 



(28) 



K{k,q,t) = K{k,q,0)e-'^^^' + 



After substituting of (j^ in (fTTjl 



P(f,t) = P^(f,t) +50 r dt'xe{rM-t')E{r,±t'), (30) 







where the upper (lower) sing corresponds to t > (t < 0) respectively. The 
total electric field is = + & ~ ~ ^ ^^"^ memory function 

Xe(r,t) = 0, t<0, (31) 

is called the electric susceptibility of the dielectric. If we are given a definite 
Xe(i) which is zero for t < 0, then we can obtain the corresponding coupling 
function /(f,a;^), 

l/(^,^fc)P = ^^^/o dtXe{r,t)smuj^t, uj^>0, 

k 

|/(f,^,^)P = 0, a;g = 0. (32) 
The operator PN{r,t) in ()3()|1 . is the noise electric polarization density 

PAr(r,t) = 

E / / d3%(a;,^,rld,(^,g-0)e-^-«*+^«^-^>/la;,^,rl4(^,g^^ 

u=i-^ w(27r)3"' 

(33) 
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Similarly after substituting in (fT5|) . we obtain the following expression 
for M(r,t), 

M(r, t) = MM{r, t) + — /'*' dt'xm{r, \t\ - t')B{f, ±t'), (34) 

/io Jo 

where Xm is the magnetic susceptibility of the dielectric 

Xm[r,t) = d^k^k\9[r,i^k)\ srncj^t, t > 0, 

Xm{r^) = 0, t<0. (35) 

If we are given a definite Xmi'l^)-, then we can obtain the corresponding cou- 
pling function g{f, lOj^) in terms of Xm{r, t) 

l5'(^,^fc)P = -r^ 2/ dtxm{r,t)smuj^t, > 0, 



A'K'^llQUj'i Jo 

\g{f,uj^,)\' = 0, u^, = 0. (36) 
The operator M^ij^t) is the noise magnetic dipole density 

M^(f,t) = 

v=i'' y(27r)3 

(37) 

Therefore the structural equations (j^ . (jHUj) and (jH^ together with the 
Maxwell equations, can be obtained directly from the Heisenberg equations 
applied to the electromagnetic field and the reservoirs. 

For a homogeneous dielectric, the coupling functions /(r, u;^) and g{f, uj^) are 
independent of r, in this case from and we deduce that Xe and Xm 
are independent of position vector r, hence the equation ()25|) can be written 

as 

V X y dt'xm{\t\ - i')V X A{r, ±t') = /^o^(r, t) + /xoV x M^(f, t), 

(38) 

where the upper (lower) sign corresponds to t > (t < 0) respectively. 
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3 The solution of Heisenberg equations 

Let us denote the fourier transform of v4(r, t) by A{q, t) 

A{f,t) = , ^ I d^qA{q,t)e"l'' (39) 

from (j2)), it is clear that 



^ [a^x{t)e^x + al^^{t)e,^x]. (40) 



The wave equation ()38p can be written in terms of A{q, t) 

i+ulA± ^ Ij*' dt'xe{\t\ - t')M ±0 - ^1 /J*' dt'xm{\t\ - t')A{q, ±t') 

= --E / ^—\u^di^^)dxik,me-'^'^%x-uJ^J*^^^^^^ 
£n N_i . /f 9'7r w 



A=l • 



where Wg- = c|g|. 

This equation can be solved using the Laplace transformation method. For 
any time dependent operator g{t), the forward and backward Laplace trans- 
formations are by definition 

gf{s) = / dtg{t)e-'\ (42) 





and 

9b{ 



s) = / dtg{-t)e-'\ (43) 
Jo 

respectively. One can solve equation (j4ip using the Laplace transformation. 
Let Xe{s) and Xm{s) be the Laplace transformations of Xe{t) and Xm(t) re- 
spectively, then Aj{q, s) and At>{.Qi s)y i-^-, the forward and backward Laplace 
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— > 

transformation of A{q,t), can be obtained in terms of Xe{s) and Xm{s) by 
taking the Laplace transformation of equation (|4ip 



£0 ''j^i^ {s±lUj^)[s^ + ul+S^Xe{s)~Ujlx^{s)] 

g%uj^)b{{k,-q,0) 



(S =F +1^1 + s'^Xeis) - UjlXm{s)] 



e^.^x} (44) 



A=l ' 



where the upper (lower) sign corresponds to Aj{q,s) {Af,{q,s)) respectively. 
Now by taking the inverse Laplace transformation of AfiQ, s) and AbiQ^ s), 

A{q, t) = r±(a;,-, t)A{q, 0) ± h±{uj^, t)Aiq, 0) 

^ 2 ^ 

— H / d^k{u^f{ujf:)^±{uj^, Uj^, t)dx{k, q, 0)e^x - ^kf*i^k)^±i^<i^ ^fc' t)d{{k, -q, 0)e^^} 
2 , ^ 

/ d-^kigiu^j^iiu^, Uj:, t)bx{k, q, 0)e^x + g*{^^i)C±{.^q, ^g, -q, 0)e^^x} 

(45) 

where 

gist 



+ \l 

^0 A=l • 



(s ± 20;^:) [s2 + a;| + s2xe(s) - u;|xm(s)] 

(46) 



where the upper (lower) sign corresponds to t > (t, 0) respectively. SOR{f{s)) 
means the inverse Laplace transformation of the complex /(s) or simply the 
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sum of residues of /(s). 

The sum of the first and second term in (j45|) is the solution of homogeneous 
equation 



A + ujlA±-j^ dt'xei\t\-t')Aiq,±t')-ulJ^ dt'xrr.i\t\-i')Aiq,±t') = Q 

(47) 

and usually tends to zero in the limit t —>■ ±00 for the case of absorbing 
dielectrics. In this case if the functions 

S + SXei-S) 



Giujg, S) 



(48) 



have poles with negative real parts, then for enough large times, 



^0 A=l 



(27r) 



d^k{ 



-dx{k,q,0) 



Ho f d^quq 



d^k{- 



'lujft+iq-T 



'UOff — U 



-6^(fc,g,0)}eVA, 



-6A(A^,<f,0) + 



(49) 



However the complete solution for A{f, t) can be obtained by substitutiing 
A{q, t) from (jUj) in ^ and using ^ 



h 



[Z±{uj^,t)e'^-'-a^x{Q) + Z*(a;,-,t)e-*«"^^4-.,(0)]e-VA 



7r)%ot^q 



^0 A=l 



d^q 



^0 A=l 



(27r)3 
^ " d^quj^ 



{2n) 



J d'HUuj^^, uo^, t)dx{k, q, 0)e^«"^^ + C±(c^,^, t)d\{k, g, 0)e-^«"^^}e>A 
d^k{r^±{uj^^, u^, t)bx{k, g, 0)e''^'' + v*^{u^„ u^, t)h\{k, q, 0)e-''^-''}e,^x, 



(50) 
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where 

Z±{uJ^,t) = SOR{^^ 2 I'^'f-^^N 2~ / \ } 

±st 

U.,..„t, = fi.,,SOR{ ^ ^^^^^^ ^ _ .^^^^^ } 

gist 

(61) 

The transverse component of the electrical field is obtained from E-^ — 
. Having the vector potential A we can easily obtain the transverse compo- 
nent of electrical polarization density and magnetic polarization density 
M using (Hp and (IMll. 

Taking the Laplace transformation of the structural equation (j3(jp . the lon- 
gitudinal component of the electrical field can be obtained as 

^11 = -— = --[ -^L= f d%Q±{u^^,t)f{u;^)ds{k,q,0)e"^'' + 

+ Q;K,,t)rK^)4(^,g,0)e-'^^1^ (52) 

where 

±st 

'''' 

In the following we consider some important examples. 
Example 1: 

Let us assume fi^j^ = g{uj%) = 0, then from and (fH^j) we have Xe(0 — 
Xmit) = 0, and from it is clear 

Z+(a;,-, t) = Z_{uj^, t) = e-'^^\ V± = C± = 0, (54) 

therefore in this case, quantization of electromagnetic field is reduced to the 
usual quantization in the vacuum, a feature that can not be seen in some of 
previous models. 
Example 2 : 

Take Xeit) and Xmit) as follows 



Xe{t) 



^ < t < A, 
otherwise, 
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otherwise 



where Xg, Xm ^ some positive constants, using 

2 _ ^c^eoXe sin^ 



(55) 

and (jHBj) we find 



k 



\9{^k\ 



hc'x 



sm^ -|- 



(56) 



and from ()30|) and ()34|) we have 



P(f,t) = P^(f,t) + 



M{r,t) = MN{r,t) + 



A J\t\-A 


Am 



t|dt'^(f,±t'), 

t\dt'B{r,±t'), (57) 



/^o(Xm + 1)A 7|t|_A 



where P/v(r, t), Mi\f{f, t) are the noise polarization densities corresponding to 
the couphng functions obtained in (fKHjl . 

In the hmit A — 0, the couphng functions tend to zero, and the relations 
(fS7|) are reduced to 



P(f,t) = £oXe^(r,t), 



M(f, f) 




Am 



/^o(Xm + 1) 



5(f,t). 



(58) 



In this limit the electrical field and the polarization densities are purely 
transverse and 



Z±{ujff, t) = cos Ugt — i 

V±i^q,t) = C±(^g,t) = 0. 



1 + Y° 

^sm u^t, 



(59) 
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The electromagnetic energy inside the dielectric is 



[^E.D + ^H.B]d'r 



^2£o(l + X°) 



+ 



2(1 + x: 



-A{q,0)-A'{q,0)] 



(60) 



where D_ is the Fourier transform of the displacement field, this energy ia 



a constant of motion contrary to the vacuum expression fl^eoE"^ + ■^]d'^i^: 



which is not a constant of motion. 
Example 3: Let Xe(i) = Pu{t) and Xm(t) = where u{t) is the step function 



u{t) 



1 t<0 
t>0 



and /? is a positive constant, then using (jH^ and (jHBj) we find 

■ 2 hc^eo/3 



l/(^. 



(61) 



(62) 



and accordingly we can rewrite ()4H) as 

A_+ujjA + PA. 



—I 



diik, -g-',0)e*-E*e%-A], 



^[d,(^,g,0)e-^-s*e>A 



(63) 



which has a dissipative term proportional to velocity. From (j3T|) one can 
obtain 

Z±{ujq;t) = e^^*[±^^ sinf^g-t + cosJ7q-t — -p:^ sin fig-t] , 



+ 



+ 



}, (64) 
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where fl^ = y^l— 4 • Therefore the stable solution of A[f,t) in the limit 
t — > ±00, is 



]e^x. (65) 



from ()52|) it is easy to show that the longitudinal component of the electrical 
field in the limit t ^ ±00 is 



_ 4(MU))e^-£*-?ng. (66) 
Example 4: A simple model for Xe{s) 

If we neglect the difference between local and macroscopic electric field for 
substances with a low density, then the classical equation of a bound atomic 
electron in an external electrical field is 

f + 7f + u'^f= -—E{t), (67) 
m 

where 7 is a damping coefficient and the force exerted on the electron due to 
atom is taken to be simply a spring force with frequency Uq and the magnetic 

force has been neglected in comparison with the electric one. If E{s) and 
r{s) are the Laplace transformations of E{t) and r(t) respectively, then from 
(jnZj) we can find 

^ --Eis) , , 

now let there be molecules per unit volume with z electrons per molecule 
such that fj electrons of any molecule have a bound frequency ojj and a 
damping coefficient 7^, then the Laplace transformation of the polarization 
density is 

hs) = ^ E 2 ks). (69) 
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2 



If ujj and 7j are identical for all of electrons, then from (j69|) 

2 _2isinz/ot , , 2 2 7 /^^^ 
Xe(i) = 2 u{t), 1^0=^0- X' 

where ^(t) is the step function defined in (jHH). We can obtain the coupling 
function /(cug) from as 

^12 _ ^c^poj^ ^ 7 7 1 

If 7 = 0, then the dielectric substance is a nondissipative one and 



l/K)l^ = ^S¥^'^(-o (72) 



In this case the noise electrical polarization density is nonzero only for fre- 
quency oj^ = uq, because Uq is the resonance frequency of equation r + UQf = 
— ■^Eqc'^'^^^ and in this frequency the energy of electromagnetic field is ab- 
sorbed by the dielectric. From (jM]) we have 



Z± a;,-, t) = SOR[ / ' 2 ] 73 



which is the sum of residues of the function inside the bracket with poles at 
±zf2-i-, where 



2 V 4 

Similarly C,± can be obtained from (|51|) with /(tu^) given in (j72j) . In this case 



(75) 

and the longitudinal component of the electrical field can be obtained from 
using f{oj^ and (5±(wg,t) given by (f72|) and (f7H|l respectively. 
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If 7 7^ 0, the substance is of dissipative kind and A{f,t), in the hmit t 
±oo can be obtained from with /(wg) given by (ffT|) and Xe(=Fza'^) 

In this case in the limit t ±oo we have 



,2 , ,2 



e--.'. (76) 



and the longitudinal component of the electrical field can be obtained from 
()52|). This example shows that this model of quantization of the electromag- 
netic field is applicable to both dissipative and nondissipative dielectrics. 



4 Concluding remarks 

By modeling the dieletric medium with two independent reservoirs, electric 
and magnetic reservoirs, we could investigate consistently the electromag- 
netic field quantization inside a linear dielectric medium. If a definite dielec- 
tric is given, i.e., Xe{t) and Xm{t) are definite functions, then we could find the 
corresponding coupling functions /(wg) and giujj:), respectively for the reser- 
voirs. In this approach, both the Maxwell and structural equations obtained 
together. In the limiting case, i.e., when there is no dielectric, the approach 
tends to the usual quantization of the electromagnetic field in vacuum, a fea- 
ture that can not be seen in some of the previous models. Also this model of 
quantization is applicable to both dissipative and nondissipative dielectrics. 
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